DO-TH-01/11 
September 2001 
revised November 20, 2001 



Nonequilibrium evolution in scalar O(N) models with 
spontaneous symmetry breaking 



Jiirgen Baacke Q and Stefan Michalski @ 
Institut fur Physik, Universitat Dortmund 
D - 44221 Dortmund, Germany 



Abstract 

We consider the out-of-equilibrium evolution of a classical condensate field and its 
quantum fluctuations for a scalar 0(N) model with spontaneously broken symmetry. 
In contrast to previous studies we do not consider the large N limit, but the case 
of finite N, including N = 1, i.e., plain A0 4 theory. The instabilities encountered in 
the one- loop approximation are prevented, as in the large- iV limit, by back reaction 
of the fluctuations on themselves, or, equivalently, by including a resummation of 
bubble diagrams. For this resummation and its renormalization we use formulations 
developed recently based on the effective action formalism of Cornwall, Jackiw and 
Tomboulis. The formulation of renormalized equations for finite N derived here rep- 
resents a useful tool for simulations with realistic models. Here we concentrate on the 
phase structure of such models. We observe the transition between the spontaneously 
broken and the symmetric phase at low and high energy densities, respectively. This 
shows that the typical structures expected in thermal equilibrium are encountered 
in nonequilibrium dynamics even at early times, i.e., before an efficient rescattering 
can lead to thermalization. 
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1 Introduction 



The investigation of the O(N) vector model at large N has a long-standing history in 
quantum field theory 0, |!|. The dynamical exploration of nonequilibrium properties 
of such models has been developed only recently [|, ^], |6], [7|, ||, [|. ^ e non ^uuibrium 
aspects of spontaneous symmetry breaking in particular have been studied in the large- N 
approximation in Refs. [HJ [TT], [12], [13], [HJ . 



At finite N there have been several approaches to formulating a Hartree-type interac- 
tion between the quantum modes. These lead, in general, to problems with renormalization 
|T5| , |T6|j . On the other hand, without such a back reaction the nonequilibrium systems with 
spontaneous symmetry breaking run into seemingly unphysical instabilities (see e.g. Fig. 9 
14]] for an illustration). For equilibrium quantum field theory at finite temperature this 



in 



interaction between quantum modes is taken into account by resummation of bubble dia- 
grams, i.e., by including daisy and super-daisy diagrams. These are essential for studying 
the phase transitions between the spontaneously broken phase and the phase with restored 
symmetry. Techniques of bubble resummation have been developed |T?|, [18|, |19|] based on 
the effective action formalism of Cornwall, Jackiw and Tomboulis (CJT) p0| . Recently 



there have been some new resummation schemes pT], p2 
renormalized. 



24| which can be consistently 



While Ref. 21 



is restricted to the lowest order two-loop graphs in the 
action, leading to bubble diagram resummation, the formalism can be extended to include 
higher order graphs p2| , ^3 , 24| . These approaches can be taken over to a formulation 
of nonequilibrium equations of motion. Here we restrict ourselves to including two-loop 
graphs of leading order in A and 1/N only; within the CJT formalism this is denoted as 
Hartree approximation [^0 . 



Having at our disposal a formalism for renormalized finite- N nonequilibrium dynamics 
we can study the new features introduced by the back reaction between quantum fluc- 
tuations in scalar quantum field theories with spontaneous symmetry breaking for finite 
N. A case of particular interest is N = 4, the 0(4) sigma model is widely used as an 
effective theory for low energy meson interactions. Its nonperturbative aspects may be an 
essential element for understanding the phenomena observed at RHIC ||25|| . Such aspects 
have been discussed previously in [^, ^7], |2S], [2D|, |3U| . Larger values of N may be realized in 
grand unified theories, whose nonequilibrium evolution may be of importance in inflation- 
ary cosmology |]3T], 



52[ |33|, |34|, |], |37|, [35], ^ ; we have included simulations for a suggestive 
value iV = 10. It should be stressed that the present approximation just constitutes a 
1/N correction, and the application to low values of N should be taken, therefore, with 



caution. Indeed, higher order corrections obtained when including the sunset diagram have 



been found to be of importance [24]] in equilibrium quantum field theory. In nonequilib- 
rium quantum field theory the role of such higher corrections is being discussed at present 

mm, Wc 



39, 40 



The out-of-equilibrium configuration that has mainly been studied is characterized 
by an initial state in which one of the components has a spatially homogenous classical 
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expectation value <fi{t). For brevity, and referring to the 0(4) case used as a model of 
low energy pion interactions, we call this component sigma (a, a = 1), the remaining 
N — 1 components pions (tt, a — 2 . . . N). In the initial state a has a classical value </>(0), 
the quantum vacuum is characterized by a Bogoliubov transformation of the Fock space 
vacuum state, a Bogoliubov transformation characterized by initial masses A^i(O) and 
A^2(0). These masses are determined self-consistently, as are their values at finite times. 

The evolution of the system is governed by the classical equation of motion for the 
field 4>(t) and by the mode equations for the quantum fields r/ a (x, t). The expectation 
values (^ a (x, t)r/ a (x, t)), a = 1. . .N, appear in both equations of motion, this constitutes 
the quantum back reaction. In the large- N limit one omits the quantum fluctuations of the 
sigma mode a = 1 and just considers the Goldstone modes. Here we are able to study the 
behavior of the a fluctuations as well, we will find that this is a rather important aspect 
in the critical region. 

The plan of this paper is as follows: in section |2| we introduce the model and the poten- 
tial obtained by bubble diagram resummation in unrenormalized form. Renormalization 
is discussed in section ||. We present our numerical results in section f| and discuss their 
implication for the phase structure of the model. Conclusions are presented in section [5]. 



2 Formulation of the model 

We consider the 0(N) vector model with spontaneous symmetry breaking as defined by 
the Lagrange density 

c= l -d^d^--^ 2 -v 2 ) 2 . (2.i) 

We consider a quantum system out of equilibrium that is characterized by a spatially 
homogenous background field. The fields are separated as 

$a = 0a(t)+7/ a (x,t) (2.2) 

into a classical part <fr a = ($ a ) and a fluctuation part r] a with (rj a ) = 0. Furthermore, in view 
of spatial translation invariance it is convenient to decompose the quantum fluctuations 
via 

ifcfct) = / {2 ^ k 2uJoi [o k /<(*,t)e <kx + 4/*(M)e-* kx ] , (2.3) 



where u 0i = yk 2 + m^. m 0i will be defined below. The subscript i = 1 denotes the sigma 
mode (a — 1), i — 2 denotes the pion modes (a — 2 ... N). 

In formulating the equations of motion and the renormalization we follow the pre- 
sentation of Nemoto et al. |21| whose generalization to the nonequilibrium system is 
st raight forward . 
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We introduce the inverse propagator in the classical background field in an O(N) sym- 
metric form 

= [ D + Mil 5 ab + ^ [M\{t) - Ml(t)} . (2.4) 

Here A4i$ are trial masses that will be determined self-consistently. In contrast to equi- 
librium quantum field theory these masses, as well as the classical field, are allowed to 
depend on time. This is not the most general parametrization of an inverse propagator; it 
is sufficiently general if the CJT formalism is restricted to the Hartree-Fock approximation, 



but not beyond it (see, e.g. [0]). If the inverse propagator has this restricted form, the 
propagator itself can be written in factorized form in terms of the mode functions fi(k,t). 
In our application the classical field has just one nonvanishing component 0i; then the 
inverse propagator has only diagonal elements and these read 

d 3 k 

■>0i 



/d k 
— — — exp («k • (x - x')) /*(*, f <) (2.5) 



(no summation over i). If the classical field has more than one nonvanishing component, 
the Green function can be likewise expressed in terms of mode functions of a coupled system 
51] , [53| . The occurrence of the mode functions allows for an interpretation in terms of Fock 



space states which has been used often to define particle numbers that refer explicitly to 
such a basis. If higher order corrections are included such naive interpretations have to be 
reconsidered. 

Superficially M.\ can be associated with the sigma mass, and M.2 with the pion mass. 
The actual meaning of these mass type parameters is more subtle, as discussed by Nemoto 
et al.. We will come back to this point later on. 



Using the ansatz ( |2.4j ) for the propagator Nemoto et al. derive the CJT effective action, 
which for a space and time independent configuration is characterized by the effective 
potential 

V^M^M,) = l -M\<f - ^ + * - 2) m 2 {Ml + (N- l)Ml} 

- 8X( N + 2) [( N + l ) M ' + 3(W - 1)^2 - 2(N - l)M\Ml + 2Nm 

+1 1 (0 ^ 2 +M2 ^^S^ ^ 2 + ^ • ^ 

Here 2 = cf) 2 . We note that our convention for the coupling constant differs from the one 
in [^TJ, furthermore we have to set m 2 = —\v 2 . 

We can easily generalize this effective potential to obtain the nonequilibrium energy 
density 

* = ^ 2 + ^0 2 -^ 4 -^^-M- M i + ( iV - 1 )- M 2} (2-7) 
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1 



(N + l)M\ + 3(JV - l)M$ - 2(N - l)M\M\ + 2NX 2 v 4 



8X(N + 2) 
+£ fl>1 + (N- l)£ fl>2 



with 



d 3 k 



\M 2 + (k 2 + Mm\fi\ 



{2.1 



(2n) 3 2uj 0i 

The equations of motion can easily be derived by requiring that this energy density be 
conserved. This is the case if the equation of motion for the field <p is given by 



4>{t) + [M\{t) - 2A0(t) 2 ] <f>(t) = , 
if the quantum modes satisfy the equations of motion 

f i (k,t)+[k 2 + M 2 i (t)]f i (k,t) = 0, 

and if the trial masses satisfy, for all t, the gap equations 

M\{t) = X [30 2 (t) -v 2 + 3F 1 (t) + (N- l)T 2 {t) 
M 2 2 {t) = X\<p 2 {t)-v 2 + T x {t) + (N + l)T 2 {t) 



Here J~i{t) are the fluctuation integrals 



d 3 k 



(2tt) 3 2^ 0? 



\fi(k,t)\'< 



(2.9) 

(2.10) 

(2.11) 
(2.12) 

(2.13) 



The gap equations incorporate the resummation of bubble diagrams. 



For a time-dependent problem we have to specify initial conditions. We choose at t — 
a value of the classical field <po = 0(0) different from its value in the equilibrium ground 
state, which is given by v apart from quantum corrections. The initial mass parameters 
m iQ = A4i(Q) are obtained by solving the gap equations (|2.11|) and (|2.12| ) at t — 0, i.e., 
by finding the extremum of the effective potential at a fixed value of <p = 0o- So the 
initial configuration is an equilibrium configuration with an externally fixed field 0o- When 
the field is allowed, for t > 0, to become an internal dynamical field the nonequilibrium 
evolution sets in. 

The equations of motion and the gap equations still contain divergent integrals and 
need to be replaced by renormalized ones. These will be derived in the next section. 

Before continuing in developing the formalism we would like to come back to the dis- 
cussion of the masses. The mass parameters M.i naively represent the effective masses for 
the a and tc fluctuations. In finite temperature quantum field theory one expects mass- 
less quanta, the Goldstone modes, if the field is in the temperature-dependent minimum 
of the effective potential in the broken symmetry phase. Likewise, in the nonequilibrium 
evolution of large-N systems it has been found [Tt], |TT|, [T2], |T^|, TA\ that the mass of the 
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fluctuations goes to zero in the broken symmetry phase as the classical field approaches 
an equilibrium value. In some sense that is trivial there, because the mass of the fluctua- 
tions appears in the classical equation of motion as well, and in equilibrium we must have 
cf) = —Ai 2 (t)4> = 0. Here the situation is more complicated and indeed we will find that 
M.\{t) will not go to zero at late times, even if the classical field settles at some constant 
value. Likewise, in the static case, at finite temperature, one finds that the "pion" mass 
M.2 is not in general zero, the question has been rised, whether Goldstone' s theorem is 
violated by the approximation or otherwise. 

The interpretation of the masses and the value of the pion mass in the present context 



have been discussed extensively in Refs.[21, 23]. These authors argue that M. 2 is not to be 



interpreted as the pion mass, but as a variational parameter which does not necessarily have 
an immediate physical interpretation. The ansatz for the inverse Green function and in 
consequence the effective potential have full O(N) symmetry. The field <fr appearing there 
is the absolute value of the field <p a . Therefore by symmetry the pion mass must trivially 
be zero in the minimum of the effective potential, which is an 0(N — 1) sphere. Nemoto et 
al. show that trivially the appropriate second derivatives of the effective potential lead to a 
vanishing pion mass. Likewise M.\ is another parameter characterizing the Green function 
and is different from the sigma mass which is given by 



ml = 2v\^-2X). (2.14) 

This is equal to 2\v 2 on the tree level only. 

It is instructive to do the simple algebra of taking second derivatives of an O(N) 
symmetric function: 

^Q- = 2cW'(0 2 ) + ^aM" (0 2 ) • (2.15) 

The pion mass is the second derivative perpendicular to the direction of 4> a , so it is given 
by the first term and vanishes if f'{<p 2 ) = 0, which defines the minimum. This corresponds 
precisely to the first derivative of the potential appearing in the classical equation of motion. 
We note that this effective mass is given by 

Ml=(M 2 l (t)-2\<p\t))- (2-16) 

It vanishes trivially if the field settles at late times at some constant value. 

The masses m 2 and M. 2 X determine the fluctuations of the classical field near the 
minimum of the effective potential in analogy to the tree level masses of the sigma and 
pion fields. If all higher corrections were included, one would expect these masses to 
determine the exact propagator near k 2 = 0. In this sense the vanishing of M. 2 \ entails a 
pole of the pion propagator at k 2 = 0. However, contrary to the large- N limit, M. 2 ^ is not 
the mass that determines the fluctuations / 2 and, thereby, the "pion" propagator Q 2 2- 
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3 Renormalization 



The renormalization of the effective potential has been discussed in [21 and [22|. As 
stated in the latter publication both formulations are equivalent, here we follow the one 
of Nemoto et al., employing the auxiliary field method, in which the counterterms in the 
effective potential are introduced via the the trial masses: 



68 = Av 2 M\ + Bv 2 M 2 2 + ~CM\ + -(N - l)DM\ . 



(3.1) 



Divergent and finite parts of the fluctuation integrals T{ and of the fluctuation energy 
densities have been analyzed in ||, using a perturbation expansion in the "potentials" 



(3.2) 



As the definition of the Green functions (|2.5| ), the expressions for the energy densities (|2.8| ) 
and for the fluctuation integrals ( |2.13| ), as well as the equations ( |2.10| ) satisfied by the 



mode functions are entirely analogous to those in ||45|| , the expansions derived there can 
be applied here in the same way. In the power series expansion of the energy densities 
with respect to Vi the zeroth, first and second order terms are UV divergent, while in the 
fluctuation integrals it is the zeroth and first order terms. In dimensional regularization 
the powers of Vi can be arranged with powers of m 2 into powers of M 2 . The expansion 
then reads 



1 



64tt 2 
1 

16tt 2 



Mf(-L + \n 



m 



iO • 



2M?m 2 



Qi 



1 4 

2 m *° 



-L + In 



m 



a, 



+ r 4 



An 



Here 



7 + In 47T 



(3.3) 
(3.4) 

(3.5) 



/i is the renormalization scale, and are the mass parameters appearing in the fluctuation 
integrals. In order to avoid bad initial singularities 57, they have to be chosen as 
m i0 = .Mj(O), the "initial masses". Finally the finite fluctuation integrals are defined by 
subtracting the UV divergent parts under the momentum integral via 



An 



d 3 k 



(2vr)32^ 0i 



\ fi \* + (k 2 + M 2 M 



(3.6) 



-2u, 



iO 



V 2 

^1 



d 3 k 



l/if-1 



Vi 
2^o 



(3.7) 



The subtractions used here are analytically equivalent to a numerically more sophisticated 
procedure used in M. The subtracted integrals are UV finite. 
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From Eqs. ( |3.3| ) and ( |3.4j ) it is evident that the divergent parts oc L are independent 
of the initial masses m 2 and thereby of the initial conditions. Exactly as in equilibrium 
theory pi 



A — B 



the divergent parts can be removed by choosing the counter terms 
and as 



as 



C = D 



1 



64vr 2 



L . 



(3.1 



With this choice the energy density as well as the gap equations become finite. The 
equations of motion are not affected. The energy density is obtained simply by replacing 
the fluctuation energies by the right hand sides of ( |3.3|) omitting the term proportional 
to L. Likewise the renormalized gap equations are obtained by replacing the fluctuation 
integrals by the right hand sides of ( |3.4j ) without the term proportional to L. In our 
numerical computations the renormalization scale \i has been taken equal to the sigma 
mass. As long as the ratio of \i and the relevant masses is far from the Landau pole, i.e., 
m 2 //i 2 <C expl67r 2 /A the dependence on the renormalization scale is weak. This is the 
condition under which the $ 4 theory can be treated as a low energy effective theory. For 
our numerical simulations with A = 1 this condition is fulfilled. 

The gap equations for the masses have to be solved once, at t — 0. By subtracting 
these equations at t — from the general gap equations one obtains the renormalized gap 
equations for the potentials Vj. These read explicitly 



V, 



A 



3(0 2 - 2 ) + 37f» + (N — l)J^ n 



16vr 2 



In 



m 



10 



Vi(*) + 



N 



fin 



+ TT + (N + 



•fin 
2 



(3.9) 



(3.10) 



16tt 2 



, ln ™?o • N+l 



] -Vi(t) + 



16tt 2 



In^o 
ft 2 



V 2 (t) 



These linear equations can be solved easily for Vi(t),i = 1,2 using a time-independent 

1 + (A/167T 2 ) ln(m 2 /mQ)) _1 appearing 



matrix. This matrix is analogous to the factor C = ( 
in the large- N case j|5| . 

The numerical implementation has been described in several previous publications (see 
e.g. 0), so we do not repeat this here. The accuracy of the computations is monitored by 
veryfying the energy conservation, which holds with a typical precision of five significant 
digits. 

We finally mention the problem of initial singularities that appears in the context of 

These can be avoided by modifying the initial quantum 



renormalization ^B), [17 
ensemble via a suitable Bogoliubov transformation. This can be done in the present model 
as well. For the values of the couplings and initial parameters considered here the initial 
singularities are numerically unimportant and have been disregarded, therefore. 
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4 Results and discussion 



4.1 Numerical simulations 

We have performed numerical simulations for the cases N = 1, N = 4 and N = 10. The 
coupling constant was taken to be A = 1 and we have varied the initial value for the field 
0o = 0(0). We have considered only values <fto > v, as for smaller values the initial mass 
17120 — A^2(0) becomes imaginary. The region <ft < v is only explored dynamically. We 
display the time evolution of the classical amplitude 4>(t) for the parameters N = 4, A = 1 
in Figs, la-c, for initial amplitudes O = 1.3v < y2v, <po = 1.445t> ~ \^2v and 0o = 1.6v > 
V2v. 




Figure la: Time evolution of the classical field cf> for O < y/2v (broken symmetry phase). Field 
amplitude and time are in units of v and v~ l , respectively. Parameters are N = 4, A = I, and 
</>o = 1.3v. 



For the tree level potential the value <fi = V2v is the value for which the energy is equal 
to the height of the barrier (N = 1) or the top of the Mexican hat (N > 1). We will discuss 
the physics associate with the three ranges of parameters below. 

In the nonequilibrium evolution of models with spontaneous symmetry breaking both 
squared masses M.\ 2 can in general take negative values. In large- N dynamics it is well 



known IT3, O, WM that in this situation the fluctations increase exponentially and drive the 



squared masses back to positive values. This stabilizes the system dynamically and prevents 
an unphysical behavior in which an exponentially increasing amount of quantum energy, 
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Figure lb: Time evolution of the classical field <p f° r 0o — v2f (critical region). Parameters 
as above but O = 1.445?;. 




200 



Figure lc: Time evolution of the classical field <\> for </> > v2v (symmetric phase). Parameters 
as above but (fi = 1.6v. 
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taken from the vacuum, is converted into classical one. Our first essential observation is 
that for all parameter sets and initial values this stabilization takes place for finite iV as 
well. We display a typical time evolution of the classical field and of the squared masses 
hA\ 2 in Fig. 2a and Fig. 2b, based on the parameters N — 4, A = 1, and <fto = lAv. In Fig. 
2a both squared masses are seen to become negative at early times and to reach positive 
values at late times. In Fig. 2b one sees the evolution of the energy of the fluctuations 
which is very strong in the time range where the squared masses are negative. 




The classical minimum of the energy is at E = — At> 4 /4, i.e., E = —.25 for this 
parameter set. So essentially all the energy is transferred to fluctuations within a few 
oscillations of the condensate field 4>. In Fig. 2c we display the two fluctuation integrals 
as functions of time for O = 2f , again for N = 4 and A = 1, i.e. in the symmetric phase. 
The pion fluctuations are seen to develop rather quickly, while the sigma fluctuations only 
develop at later times. 

The main aspect we will consider here is the phase structure of the model. In ther- 
mal quantum field theory one expects a phase with spontaneously broken symmetry at 
low temperature, and symmetry restoration at high temperature. Here we have a micro- 
canonical description, so instead of temperature we specify the energy which in turn is 
determined by the initial value 0o- As already mentioned the system will not thermalize 
in the approximation used here, but it is characterized at late times by limiting values, or 
time averages, attained by the various masses, and by the classical field <fi. The latter can 
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Figure 2b: Time evolution of the energy densities. Parameters as in Fig. 2a. Displayed are 
the classical energy density, the fluctuation energy density including finite renormalization terms 
and the total energy density. 




Figure 2c: Time evolution of the fluctuations integrals. Parameters as in Fig. 2a, however with 
4>q = 2v. upper curve: pion fluctuations, lower curve: sigma fluctuations. 
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be compared with the temperature-dependent vacuum expectation value v(T), the former 
ones are associated to correlation lengths. 




Figure 3a: The late-time amplitude as a function of the initial amplitude. Parameters are 
N = 1 and A = 1. 



Classically the spontaneously broken phase is expected to be situated in the interval 
v < 0o < V2v. For larger values of O one should be in the phase of restored symmetry. 
Since the effective potential receives quantum corrections the boundaries of these intervals 
will be slightly shifted, for brevity we continue to use the classical values in the following 
discussion. 



4.2 Late time behavior of the classical field 

We first consider 0oo(0o), the value of averaged at late times, which can be considered as 
an order parameter, replacing the vacuum expectation value v(T) of finite temperature field 
theory. In the analysis of the broken symmetry phase in the large- N limit it was observed 
|12| , p~3| , jllf that the value of at late times has a specific and universal dependence on the 
initial value 0o, given at zero temperature by 



(4-1) 



We find that this is the case for finite iV as well. The dependence 0oo(0o) is shown in 
Figs. 3a-c for N — 1,4 and 10, respectively, the curves are rough fits of the form of 
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Eq. |4.1| , taking into account a slight shift of the point of "phase transition" ~ \f2v. 
Such a functional behavior is typical of a phase transition of second order. In a first oder 
phase transition the vacuum expectation value of the true minimum jumps at the phase 
transition from the broken symmetry minimum to the symmetric one. As the functional 
behavior of Eq. ( |4.1| ) has an infinite slope at the critical point it is hard to decide from 
numerical data whether 0^ near 0q = 

y/2v goes to zero continously or via a discontinuity. 
The results strongly indicate that the transition is discontinuous. This is suggested as well 
by the time evolution of displayed in Fig. lb. There seems to be a minimum at = 
and another one near ~ .5v , as one knows it from typical finite temperature potentials 
of first order phase transitions. However, even if the transition is first order, it is very close 
to a second order one. In thermal equilibrium a first order phase transition is expected 
2M, it becomes second order only after including higher loop corrections. Fig. lb 



puts in evidence that the time averaging is problematic near the phase transition; this 
implies that the functional dependence the "order parameter" 0^ near 0o ~ \/2v is not 
determined with high precision. 

We would like to remark that we have chosen special initial conditions which allow the 
classical field to move in one fixed direction only. So the system cannot see the differ- 
ence between a double-well and a Mexican hat potential. Imposing more general initial 
conditions |53] may lead to an improved understanding of the nonequilibrium properties 



of the system. Such a generalization may be useful especially near the transition point 
phi ~ \/2v. 

For > \plv the time average of vanishes as t — > oo, and in fact already at an early 
stage of evolution. However, the amplitude of the oscillations decreases very slowly, if at 
all. So, although the order parameter 0^ shows the behavior expected in the symmetric 
phase of a thermal system, and although the system becomes essentially stationary, it 
shows no resemblance to a system in thermal equilibrium with a constant value of 0. We 
will further analyze this phase, below. 



4.3 Late time behavior of the masses 

Another set of variables characteristic for the phase structure of the O(N) model in equi- 
librium are the mass scales or correlation lengths. In the broken symmetry phase we expect 
a vanishing pion mass and a nonzero value of the field 0. In section [| we have discussed 
to some extent the problem of the Goldstone mass, which should not be identified naively 
with M.2-, rather, we have proposed to associate it with the "classical mass" M c \ defined 
in Eq. ( |2.16|) . Both masses would coincide in the large- N limit. 



Likewise, the sigma mass is not given naively by Mi- We plot the three squared masses 
M\, Ml, and Mix as functions of O in Fig. 4, for N = 4, A = 1. M\ is seen to be almost 
zero for v < 0o < V%v and to increase for larger values of O , as expected for the Goldstone 
mass. The mass M\ is small for v < 0o < \/2v but definitely different from zero, and 
increases likewise for v > y/2v. The "sigma mass" Mi is different from zero everywhere, 



14 




1,2 1,4 1,6 1,8 



Figure 4: The various masses averaged at late times (t=300), as functions of the initial 
amplitude. Squares: M.\, diamonds: M\, circles: M.1 V Parameters are N = 4 and A = 1. 



except near 



V2z 



The formalism presented here contains of course the limit N — > oo. This limit is 
obtained by letting A oc 1/N and 0, 0o, v oc y/N. In this limit the quantum fluctuations r\\ 
become irrelevant, and so does M.\. The "pion" mass M.2 and the "classical mass" M. c \ 
become identical to the mass M. of Refs. \T1, |T^, |TB|, 14] or x(t) of Ref. 



Toil . In the broken 



symmetry phase the instability pushes this mass to zero. Here, for finite N, we find Ai 2 . to 
remain positive though small in the broken symmetry phase, while Ai^i again vanishes. It 
should be mentioned that the vanishing of the classical mass M.\ in the broken symmetry 
phase is due to the fact that the field <fi(t) settles at a constant nonzero value. For the 
fluctuation mass M.\ the dynamics of backreaction only forces it to be nonnegative. In the 
large-iV limit the identity of the two masses Ai 2 = A\\ = M. 2 X entails their vanishing at 
late times, as expected from the Nambu-Goldstone theorem. The situation found here at 
large times is analogous to the one found in thermal equilibrium |21, |23 



(see section 



For a phase transition of second order the a mass, all mass scales are expected to vanish 
at the "critical point", i.e., for <pQ ~ \[2v. We see from Fig. 4 that this is almost the case. 
For a first order phase transition the curvature of the potential in the broken symmetry 
minimum decreases when approaching the phase transition, but remains positive up to and 
beyond the phase transition. Above the phase transition the curvature in the symmetric 
minimum increases again. The fact that M\ does not really reach zero may of course be 
a deficiency of our numerics. We point out, however, that critical behavior implies large 
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length and also time scales, so that in this neighborhood the time averaging is precarious 
and the late time values are not very precise. It is not clear whether the field in Fig. 
lb will jump back and forth between the various minima at later times again and, if it 
does, whether this is not a consequence of accumulated tiny numerical errors. We again 
conclude that, if the phase transition is of first order, it is "weakly first order". While we 
have mostly presented results for N = 4, we find similar results for N — 1 and N — 10. 
For N — 1 the mass M.2 is of course meaningless. 

As we have mentioned above, in the symmetric phase the time average of the field (f) 
becomes zero at late times, while the field itself continues to oscillate with essentially 
constant amplitude, implying that the time average (<f) 2 ) of the squared field remains 
different from zero. This implies that the masses M.\ and Ai 2 necessarily have different 
time averages; indeed the classical field contributes with 3A(0 2 ) to (Aif) but only with 
A(0 2 ) to (A4 2 ,). Actually even the quantum fluctuations T\ and Ti have different time 
averages. This is contrary to what one expects in a symmetric equilibrium phase. The 
problem is inherent in the initial conditions which necessarily require a high excitation and 
the choice of direction in O(N) space. It is appropriate at this point to compare with a 
$ 4 model without spontaneous symmetry breaking, i.e., with m 2 = — Xv 2 > 0. There again 
the classical field oscillates around = with slowly decreasing amplitude. 
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Figure 5: The various masses averaged at late times (t=300), as functions of the initial 
amplitude for a model without symmetry breaking, N = 4, A = 1, v 2 = — m 2 /A = — 1. 
Symbols as in Fig. 4. 

Again the time average of mass M.\ remains different from the one of Ai 2 pertaining to 
the fields rj a with a = 2 . . . N. We show in Fig. 5 the dependence of the late time averaged 
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masses for this manifestly symmetric theory. The close resemblance with the behavior in 
the high energy phase of the model with spontaneous symmetry breaking is obvious. The 
difference between the masses is not related to spontaneous symmetry breaking but to the 
fact that the system is prepared in a non-symmetric state. 

While the behavior of masses J\A\ and M. 2 C \ found here essentially reproduces the one 
of J\4 2 in the large- N limit we have seen here, that the sigma mass contains valuable 
information about the nature of the phase transition, not available in the large- N limit, 
and that the sigma fluctuations are relevant for finite N. 



4.4 Correlations 

Correlations of mode functions have been discussed in various publications [|H], |28], [13|, ^9] . 
If one thinks of the N = A model as a model for pion production, possibly displaying struc- 
tures like disordered chiral condensates, it is the correlations between the pion fluctuations 
{i = 2 or a = 2 . . . N) which are relevant. In the large- N limit these are the only available 
correlations. It was found in a large-iV computation including the full backreaction 



that the correlation length grows with time in the broken symmetry phase. Kaiser p8 



analyzing the correlations in the initial phase, before back reaction sets in, suggested the 
occurrence of large correlation lengths. Hiro-Oka and Minakata |29j performing a computa- 



tion with backreaction, however unrenormalized and with a different Hartree factorization, 
suggested that the correlation length remains small. 

The pion correlations are obtained as the Fourier-Bessel transform 

/cl k 
, 9 ^3 9 ^\f 2 (k,t)\ 2 

= — / dk ksin(kr)\f 2 (k,t)\ 2 . (4.2) 
2ir z r Jo 

We have performed a series of simulations illustrating the transition from N — > oo to finite 
N. Fig. 6a shows the correlation functions rC(r) at times t = 30, 50 and 100 for X-N = 1, 
with iV = oo, 10 and 4. 

The initial amplitude is 4>q ~ 1.024t> implying in all cases m 2 Q — 0.2. It is seen that 
there are long range correlations in all cases. They are seemingly related to or actually 
generated by the oscillations of the classical field. Its amplitude is displayed in the same 
figure for the case N — > oo, the finite iV amplitudes behave very similarly. The correlations 
are positive throughout for the large-iV case. For iV = 10 and N = 4 they alternate in 
sign in the central region at late times. So at small N, in particular at iV = 4, these 
results certainly do not suggest a growth of positively correlated domains but rather a 
wash-out implied by the alternating signs. In fact the size of a positively correlated region 
is essentially of the order of an oscillation period of the sigma field. 

The close relationship between the oscillations of the classical field and the correlations 
is further demonstrated in Fig. 6b. For N = 4, X = 1 and an initial value 0o — l-2i> 
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Figure 6a: Equal-time correlations of pion fluctuations. We plot rC{r), with C{r) as defined 
in Eq. flO]), at times t = 30, 50, 100, for parameters v = VN, X = 1/N, N = 4, 10 and for 
N — > oo. Initial amplitude O — l-024t>. On top the amplitude <f>(t)/v for N — > oo. 
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Figure 6b: Equal-time correlations of pion fluctuations; parameters N = 4, A = 1/4, v = 2, 
initial amplitude (fi = 1.2v. 
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the oscillations of the classical field decay rather quickly, and so do the correlations. This 
means that one has to be cautious with their interpretation: these correlations are not 
due to an interaction between field fluctuations but are generated by the coherence of 
such fluctuations with an external source. This is not necessarily the wrong physics, but 
certainly this feature is inherent in the mean field approach. 



4.5 Momentum spectra 

In the large- N limit one of the most pronounced characteristic of the momentum spectra 
of the "pion" fluctuations is the occurrence of parametric resonance bands. These de- 
velop already in the early stage of evolution, before backreaction sets in. Then the time 
dependence of the classical field is described by Jacobian elliptic functions, and the fluctu- 
ations are solutions of the Lame equation. These solutions have been derived and discussed 
extensively in [|TT], [TJ, [T3[ . 



For our finite-iV system the early time behavior of the classical field and of the pion 
modes is the same as in the large- N limit, as before the onset of fluctuations both Ad^i 
and M.2 are approximatley equal to A(0 2 — v 2 ), while M.\ ~ A(30 2 — v 2 ). The analysis 
of the pion fluctuations can be taken over from Ref. [11]], therefore. For O < V%v the 
parametric resonance band is given by [] 




(4.3) 

as given in Ref. [IT]. F° r 0o > parametric resonance occurs (see Appendix A) in the 



momentum interval 

< k 2 < . (4.4) 
In Figs. 7 and 8 we plot the integrands of the various fluctuation integrals J-'i, i.e. 

^M [k - t)? - l+ £i) (45) 

as functions of k. 

As seen in Figs. 7b and 8b parametric resonance indeed dominates the pion fluctua- 
tions within the momentum intervals ( |4.3| ) and ( |4.4| ) in the broken symmetry and in the 
symmetric phase respectively. The sigma fluctuations are unimportant in the broken sym- 
metry phase, as displayed in Fig. 7a. In the symmetric phase they develop a pronounced 
resonance which seemingly becomes sharper and sharper as time increases. The position of 
the resonance can be obtained by considering the periodicity T of the square of the classical 



3 Here and in the following we simplify the discussion by neglecting modifications due to finite renor- 
malization corrections. 
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Figure 7a: Spectrum of sigma fluctuations in the broken symmetry phase. We display the 
integrand defined in Eq. ( |4.5| ) of the fluctuation integral as a function of k for A = 1, v — 1, 

0o = 1.2u at t = 100. 
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Figure 7b: Spectrum of pion fluctuations in the broken symmetry phase. Notation and param- 
eters as in Fig. 7a. 
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Figure 8a: Spectrum of sigma fluctuations in the symmetric phase. Notation and parameters 
as in Fig. 7a, however O = 2v. 
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Figure 8b: Spectrum of pion fluctuations in the symmetric phase. Notation and parameters 
as in Fig. 7a, however O = 2v. 
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field and the frequency of the "sigma" fluctuations; one expects k 2 + (A4 2 ) — 47r 2 /T 2 . The 
empirical ("measured") values of T and (Aif) for the simulation with our parameter set 
yield k 1.5 in agreement with Fig. 8a. As the resonance is due to a time-dependent mass 
term, it should be situated in a resonance band, as indeed suggested by Fig 8a. For this 
resonance band the early time analysis does not apply, as the sigma fluctuations evolve at a 
time when backreaction due to pion fluctuations has already modified the equation for the 
condensate field </>(£), see Fig. 2c. It should be noted that an analoguos resonance appears 
in the spectrum of sigma (longitudinal) fluctuations in the <3> 4 model without spontaneous 
symmetry breaking, while the pion (transverse) fluctuations display a broad parametric 
resonance band. So again the difference between the fluctuations parallel and orthogonal 
to the classical field is not related to spontaneous breaking but to the non-symmetry of 
the state. 



4.6 A naive dynamical effective potential 

As we have seen, the behavior of the system below 0o = V%v displays the features of a 
spontaneously broken phase in a rather convincing way, while for the symmetric phase the 
non-symmetric initial conditions remain manifest even at late times. The amplitude of 
oscillations of the classical field does not go to zero even if we extend the simulation to 
times beyond which we believe our numerical reliability. 




The evidence for symmetry restoration remains indirect, therefore. So it might be 
useful to have yet another criterion. This is given by the dynamics of the classical field. 
In a rather trivial way the potential that yields the "force" experienced by the classical 
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Figure 9b: Evolution of the potential energy. Parameters as in Fig. 9a, however O = 2t>. 



field is given by the potential energy, obtained by subtracting the kinetic energy from the 
total energy, V pot = E — <p 2 /2. This potential is of course time-dependent and is explored 
only within the range of oscillations of the field. We display this potential energy in Figs. 
9a and 9b for the broken and for the symmetric phase, respectively. It is seen that the 
positions of minima shift from = ±v to smaller absolute values in the broken symmetry 
phase, while in the symmetric phase the separate minima entirely disappear after a few 
oscillations and a new minimum appears at = 0. 




5 Conclusions 

We have presented here an analysis of nonequilibrium dynamics in O(N) models with 
spontaneous symmetry breaking at finite N in a bubble-resummed one-loop approximation. 
This has allowed us to study some new features of such systems, not accessible to the large- 
N or one-loop approximations. 

We have found that the back-reaction of the quantum fluctuations on themselves pre- 
vents, as in the large- iV approximation, a catastrophic instability of the system as found in 
the one-loop approximation. The mechanism is, like in large-iV, an exponential evolution 
of those quantum fluctuations whose effective mass squared has become negative, pushing 
this mass squared back to positive values. 

The dependence on the initial conditions displays features of an expected phase tran- 
sition between a regime with spontaneous symmetry breaking and a symmetric phase. A 
more detailed analysis of nonequilibrium dynamics in the critical region may provide new 
insights into the nature of the phase transition. 
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Our approach provides an new tool for performing studies of nonequilibrium evolution 
in scenarios typical of heavy ion collisions |25j or of the dynamics of inflation in the early 
using appropriate models like the 0(4) sigma model or realistic grand unified 



universe 
theories. 
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It should be especially useful for studying phase transitions and critical behavior 
in an nonequilibrium context. It would be useful to extend this study to more general 
initial conditions exploiting more of the mult i- dimensional mexican hat structure of the 
potentials in models with N ^ 1. The techniques using coupled-channel Green functions 



are available 



Thermalization is not expected in the large- N approximation and, in the absence of a 
rescattering involving the next order sunset diagrams, may not be expected here either. 
One may consider this as a drawback of the Hartree type approximation we are using. In 
simpler quantum systems where one can compare to the exact time evolution the Hartree 
approximation is found to constitute an improvement with respect the large- N approx- 
imation, but is further improved by the inclusion of sunset diagrams P5| , |55| , pEOfl . Still 
it describes well the early time behavior. In heavy ion collisions or in the early universe 
the behavior at early times, and/or without reaching equilibrium and thermalization, may 
be even more realistic and is therefore interesting as well. It is rewarding that the phase 
structure of the theory is revealed by nonequilibrium dynamics at early times already. 

Thermalization is, on the other hand, a basic theoretical issue, independent of realiza- 
tion in concrete physical processes. It has recently been investigated by studying the time 
evolution of classical fluctuations in Hartree approximation |^T[ . As the present formalism 
can be extended to include higher-loop diagrams |22, ^3], |24| our work can be considered as 
a first step towards calculations including rescattering of fluctuations [32 and of "controlled 
nonperturbative dynamics" |43, [44j] in three dimensions, using a continuum regularization. 
This requires the inclusion of sunset and higher order diagrams, incorporating rescattering 
of the quantum fluctuations. Such calculations are being performed at present p2| , using 
lattice regularizations. We think that it will still take a long time until the limitations 
introduced by various approximations are fully understood, leading to a well-based under- 
standing of thermalization in quantum field theory. So various alternative approaches will 
have to be considered. With our modest new step, introducing some 1/N corrections we 
are clearly still far off such a demanding formal and in particular numerical task; we think, 
however, that our investigation provides some useful, and possibly inspiring, new insights. 
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A Some solutions of the Lame equation 



The following analysis follows closely the one in Ref. [p|. The basic formulae are from 
Ref. pi] , the only difference in notation is in the argument of the elliptic integral K which 
is m there and k = y/m here. We introduce the dimensionless classical field x — 4>l v an d 
the dimensionless time variable r = y/Xvt. We denote differentiation with respect to r by 
a prime, df jdr = /'. Then the equation of motion for x before the onset of backreaction 
via fluctuations, and neglecting finite renormalization terms | is given by 



X" + x 3 



X 



. 



(A.l) 



Using the dimensionless momentum variable q = kj y/\v (here k is the momentum variable 
used in the main text) the mode equations before the onset of backreaction read 



/r + (g 2 -l + 3 X 2 )/i = 

fi + (q 2 -i + x 2 )f2 = o. 



(A.2) 
(A.3) 



For the symmetric phase, i.e., Xo — 4>o/ v > v2, the solution of the classical equation of 
motion is given by 

X(r) 



Xo 



1 -sn 2 (rJx 2 o 



(r y^l, k) 



1/2 



(A.4) 



where sn denotes the Jacobi elliptic function, its index k is given by k = xo/ \/2(Xo — !)• 
It may be related to a WeierstraB elliptic function p via 



sn 2 0Vxo - 1, k) 



1 



k 2 (el - e3) 



P 



X 2 -l + iK'(k) 



- e 3 



(A.5) 



Furthermore we have (see section 18.9 of Ref. iK'(k) = o/i/ei — e3 where u;' is the 

imaginary half period of the WeierstraB function p. For the roots we find e\ = Xo/2— 2/3, 
62 = 1/3, and e-s = —Xo/^ + 1/3; the invariants are given by g 2 = (xl — l) 2 + 1/3 and 
3g 3 = — (Xo — l) 2 + 1/9. The half periods of the double-periodic function p are related to 



the roots by p{oS) = ex, p(oo + u') = e 2 and p(u') 
then become 



e$. The equations of motion for the fi 



/r+[g 2 + l-6p(r + ^)]/i 
1 



fS + 



2p(r + cu') 



h 





. 



(A.6) 
(A.7) 



The general Lame equation reads 55 



/"-[a + 7i(7i + l)p]/ = 0. 



{hi 



4 This simplification is not really necessary, but appropriate for the couplings used here. 
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For f 2 we have n = 1 and the solution is given by 



/2= o(r±a) exp 
o-(r) 

where a is determined by the equation 

p(a) = a , (A. 10) 

and where a and ( are WeierstraB functions. 

For fi we have n = 2 and the solution is given by []5^] 

= d a(r± a ) exp(Tr(c(a) + ^ (A n) 

ar (TIT 



where a is one of the solutions of 



with ai = a/3, and where 



P(a) = (A.12) 
2>a{ - g 2 



P= . (A.13) 

2p(aJ — ai 



The solutions of the Lame equation are quasiperiodic; if the argument increases by a 
period 2a; then the solution reproduces itself up to factor exp(iF(a)) where F(a) is known 
as Floquet index. If it is imaginary then there is an exponentially increasing solution, if it 
is real the solutions are periodic up to a phase. The right hand sides of ( |A.12| ) and ( |A.10| ) 



are real. The WeierstraB function p maps the fundamental rectangle [0, u>, uo + uo' , oo', 0] onto 
the upper half plane, so the solutions a of these equations are situated on the boundary of 
the fundamental rectangle; the origin a = is mapped to the infinite point. 

For f 2 a closer analysis shows that on the sections [0, u\ and [u/, u' + u] of this boundary 
the Floquet index 

F(a) = -2i [aCM - wC(«)] ( A - 14 ) 

is imaginary, so in these regions one has parametric resonance ( "forbidden bands" in anal- 
ogy to solutions of the Schrodinger equation in periodic potentials); on the sections [0, a/] 
and [uj,u> + uj'] the solutions are oscillatory ("allowed bands"). The "forbidden" bands are 
characterized by t\ < p(a) < oo and < p(a) < e 2 . For f 2 this implies parametric reso- 
nance in the in the intervals — oo<g 2 <— %q/2 + 1<0 which is excluded kinematically, 
and in the interval < q 2 < This resonance band manifests itself in Fig. 86. 

For fi, the sigma fluctuations, the analysis is somewhat more cumbersome, Eq. ( |A.12j ) 
becomes explicitly 

P(a) = "9 (1 +9 T-1-3(X?-D 2 ■ <A - 15) 
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Again a is situated on the boundary of the fundamental rectangle with the same corners. 
The Floquet index is found to be 



^'^'-^-" wfa ' (A - 16) 

It is imaginary for values of a on the sections [0, u] and [u', uj' + u] of the boundary of 
the fundamental rectangle. Again this maps onto the intervals e\ < p{ot) < oo and e 3 < 
p(a) < €2- The analysis of the function on the right hand side of Eq. ( |A.15| ) shows that 



the resonance bands are — oo < q 2 < — 1 — j3(xo — l) 2 + 1 and — 3xo/2 < q 2 < which are 



excluded kinematically, and, in the physical region, 3xo/2 — 3 < q 2 < — 1 + y 3(xo — l) 2 + 1. 
The latter resonance band is not manifest in Fig. 8a as the sigma fluctuations develop only 
after the classical equation of motion is modified by the backreaction due to the pion 
fluctuations. We have verified that our analytical result for the resonance band is correct 
by running the simulation without backreaction. The parametric resonance of the sigma 
fluctuations analyzed here may manifest itself for other parameter sets and the result may 
be of importance, therefore. 
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